Random Delays and the Synchronization of Chaotic Maps 
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We investigate the dynamics of an array of logistic maps coupled with random delay times. 
We report that for adequate coupling strength the array is able to synchronize, in spite of the 
random delays. Specifically, we find that the synchronized state is a homogeneous steady-state, 
where the chaotic dynamics of the individual maps is suppressed. This differs drastically from 
the synchronization with instantaneous and fixed-delay coupling, as in those cases the dynamics is 
chaotic. Also in contrast with the instantaneous and fixed-delay cases, the synchronization does not 
dependent on the connection topology, depends only on the average number of links per node. We 
find a scaling law that relates the distance to synchronization with the randomness of the delays. We 
also carry out a statistical linear stability analysis that confirms the numerical results and provides 
a better understanding of the nontrivial roles of random delayed interactions. 
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A system composed of many nonlinear interacting 
units often forms a complex system with new emergent 
properties that are not held by the individual units. Such 
systems describe a wide variety of phenomena in biology, 
physics, and chemistry. The emergent property is usu- 
ally synchronous oscillations. In biology examples are the 
synchronized activity in pacemaker heart cells, the cicar- 
dian rhythms that control 24h daily activity, and the 
flashing on-and-off in unison of populations of fireflies. 
Nonbiological examples include synchronized oscillations 
in arrays of electrochemical reactions, in laser arrays, and 
in Josephson junction arrays fl']. 

The effect of time-delayed interactions, which arise 
from a realistic consideration of finite communication 
times, is a key issue that has received considerable atten- 
tion. The first systematic investigation of time-delayed 
coupling was done by Schuster and Wagner Q , who stud- 
ied two coupled phase oscillators and found multistabil- 
ity of synchronized solutions. Since then, delayed inter- 
actions have been studied in the context of linear sys- 
tems , phase oscillators , limit-cycle oscillators J3 , 
coupled maps , neuronal 0, 0, lioj and laser ^ 
systems. Most studies have assumed that all the inter- 
actions occur with the same delay time (only few have 
consider non-uniform delays 

113, 111 El)- However, ac- 
tual delays in real extended systems are not necessarily 
the same for all the elements of the system, they might 
be distant-dependent or randomly distributed. In popu- 
lations of spatially separated neurons, the synaptic com- 
munications between them, which depend on the prop- 
agation of action potentials over appreciable distances, 
involve distributed delays. In computer networks, ran- 
dom delays arise from queueing times and propagation 
times. In epidemic dynamics, migrations of geographi- 
cally spread populations lead to different delays which 
infiuence the transmission of diseases. 

While it is well-known that oscillators which interact 



with different delay times can synchronize (an example is 
the synchrony arising in different neuronal groups of the 
brain which might lead to both, epilepsy and Parkinson 
disease), the mechanism by which this synchrony arises 
and the infiuence of the random communication times 
remains poorly understood. The focus of this letter is 
to investigate the infiuence of such random delays in the 
synchronization of a simple model of coupled chaotic os- 
cillators. We consider an ensemble of logistic maps and 
show that, in spite of the random delays, for adequate 
coupling strength the array is able to synchronize. Sur- 
prisingly, in the synchronized state the chaotic dynamics 
of the individual maps is suppressed: the maps are in a 
steady-state, which is unstable for the uncoupled maps. 
This is in sharp contrast with the cases of instantaneous 
and fixed-delay coupling, as in those cases the dynamics 
of the array is chaotic. By studying the transition from 
chaotic synchronization to steady-state synchronization 
as the randomness of the delays increases we discover a 
scaling law that relates the distance to synchronization 
with the randomness of the delays. We also investigate 
the infiuence of the array topology and find that steady- 
state synchronization depends on the average number of 
links per node but not on the array architecture. This is 
also in contrast with the instantaneous and fixed-delay 
cases, as in those cases the synchronization depends also 
on the connection topology 0- Finally, we present a 
statistical linear stability analysis that demonstrates the 
stability of the solution found numerically. 

We consider the following ensemble of N coupled maps: 



N 

+ 1) - (1 - e)/[x.(i)] + ^ ^ ri,,f[xj{t - r,,)] (1) 

Here i is a discrete time index, z is a discrete spatial 
index (i = 1 . . . N), f{x) = ax{l — x) is the logistic map. 
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FIG. 1: Xi vs. e (a) and (b); vs. e, (c) and (d). In Figs. 
1(a) and 1(c) the delays are distributed exponentially (tq — 0, 
c = 1.2, see text); the distribution is shown in Fig. 1(e). 
In Figs. 1(b) and 1(d) the delays are Gaussian distributed 
{to = 3, c = 2, see text); the distribution is shown in Fig. 
1(f). The inset in Fig. 1(d) shows with detail the transition 
to synchronization; decreases abruptly at e ~ 0.4, and 
is zero for e > 0.8. Parameters are: A'^ = 500, a = 4, and 
p = 0.3. 



the matrix 77 — {rjij) defines the connectivity of the array: 
Vij — Vji = 1 if there is a hnk between the ith and jth 
nodes, and zero otherwise, e is the couphng strength 
and Tij is the delay time in the interaction between the 
ith and jth nodes (the delay times and Tji need not 
be equal). The sum in Eq.JQ) runs over the bi nodes 
which are coupled to the ith node {hi — 77^ ) . The 
normalized pre-factor l/bi means that each map receives 
the same total input from its neighbours. 

First we note that the homogeneous steady-state 
Xi{t) — Xj{t) — xq, V i, j, t, where xq is a fixed point of 
the uncoupled map, xq — /{xq), is a solution of Eq.QJ 
regardless of the delays and of the connectivity of the 
array. 

Next, let us present some results of simulations that 
show that this state, with xq being the nontrivial fixed 
point, Xq = 1 — 1/a, can be a stable solution for adequate 
coupling and random enough delays. The simulations 
were done choosing an initial configuration, Xi{0) ran- 
dom in [0,1], and letting the array evolve initially with- 
out coupling [in the first time interval < i < max(Ty )]. 
We present results for = 4, corresponding to fully de- 
veloped chaos of the individual maps, but we have found 
similar results for other values of a. We illustrate our 
findings using the small- world topology ) but we have 
found similar results for random and regular topologies 
fl5j . as discussed below. 

With both, either randomly distributed or fixed delay 
times, if the coupling is large enough the array synchro- 
nizes in a spatially homogeneous state: Xi{t) = Xj{t) V 
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2: Xi vs. e (a) and (b): 
:0 V j,i; in (b),(d) r,, = 



a'^ vs. e (c) and (d). In (a),(c) 
3 V All other parameters as 



in Fig. 1. 



i, j. Figures 1 and 2 display the transition to synchro- 
nization as e increases. At each value of e, 100 iterates of 
an element of the array are plotted after transients. To 
do the bifurcation diagrams we varied only e; the connec- 
tivity of the array {r]ij), the delays {Tij), and the initial 
configuration [a:i(0)] are the same for all values of e. Fig- 
ure 1 (a) displays results for delays that are exponentially 
distributed. Fig. 1(b) for Gaussian distributed, Fig. 2(a) 
for zero delay, and Fig. 2(b) for constant delays. It can 
be observed that for small e the four bifurcation diagrams 
are similar; however, for large e they differ drastically: Xi 
is constant in Figs. 1(a) and 1(b), Xi = xq — 1 — 1/a, 
while Xi varies within [0,1] in Figs. 2(a) and 2(b). 

To characterize the transition to synchronization we 
use the indicator tr^ = 1/7V < ^i[xi{t)— < x >]^ >t, 
where < . > denotes an average over the elements of the 
array and < . >t denotes an average over time. Figures 
1(c), 1(d), 2(c) and 2(d) display vs. e for the bifurca- 
tion diagrams discussed above. It can be observed that 
for large e there is inphase synchronization in the four 
cases {xi{t) = Xj{t) and ct^ — 0]; however, we remark 
that an inspection of the time-dependent dynamics re- 
veals that for randomly distributed delays the maps are 
in a steady-state, while for fixed delays the maps evolve 
chaotically. It can also be observed that the four plots 
vs. e are similar for small e (in Fig. 2(d) the array 
synchronizes also in a window of small e; this occurs for 
odd delays and was reported in 0)- 

Let us now investigate the transition from chaotic syn- 
chronization (for fixed delays) to steady-state synchro- 
nization (for random delays) by introducing a disorder 
parameter c that allows varying the delays from constant 
to distributed values. Specifically we consider 

i) '''ij = ''0 + Near(c^), where ^ is Gaussian distributed 
with zero mean and standard deviation one and Near 
denotes the nearest integer. The delays are constant 
{Tij = To) for c = and are Gaussian distributed around 
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FIG. 3: (a) vs. c and (b) vs. c*. The delays are 
Gaussian distributed with ro = 2 (o), 3 (A), 4 (V), 5 (<) ; 
exponentially distributed with ro = 3 (x), 4 (+), 5 (*), 6 (o). 
e = 1, all other parameters are as in Fig. 1. 

To for c ^ [depending on tq and c the distribution of 
delays has to be truncated to avoid negative delays, see 
Fig. 1(f)]. 

ii) Tij = To + Int(c^), where ^ is exponentially dis- 
tributed, positive, with unit mean and Int denotes inte- 
ger. The delays are constant (t^ = tq) for c — Q and are 
exponentially distributed, decaying from tq for c 7^ 0. 

Figure 3(a) displays the transition between the two 
synchronization regimes as the randomness of the delays 
increases. We plot cr^ vs. c for different delay time dis- 
tributions and e large enough such that the array syn- 
chronizes for c = 0. It can be observed that the array 
desynchronizes as c increases and the delays become dif- 
ferent from each other. There is a range of values of 
c such that the delays are not random enough to induce 
steady-state synchronization; however, for c large enough 
the array synchronizes again and tr^ = 0. 

To investigate this transition we considered a normal- 
ized disorder parameter c* — Dr/ < t >, where < r > 
is the average delay and Dt is the standard deviation of 
the delay distribution. By plotting ct^ vs. c* [Figs. 3(b)] 
we uncover a scaling law: as < r > increases the curves 
collapse into curves of similar shape, and the transition 
to steady-state synchronization occurs for c* ^ 0.5. 

The value of the disorder parameter above which the 
arrays synchronizes in the steady-state depends on the 
coupling strength. Figure 4(a) displays the synchroniza- 
tion region in the parameter space (e, c). To determine 
the synchronization region we did simulations with differ- 
ent initial conditions, array connectivities, and delay time 
realizations: the black indicates parameters for which 
the array synchronized in all the simulations, while the 
white indicates parameters for which the array did not 
synchronize in any of the simulations. The grey region 
in the boundary of the synchronization region indicates 
that there are some initial conditions and/or realizations 
of rjij and Tij for which the array did not synchronize. 
Two different synchronization regions can be clearly dis- 
tinguished: for c — Q and for c > 0.5. The former cor- 
responds to chaotic synchronization for fixed delays, and 
the latter, to steady-state synchronization for distributed 
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FIG. 4: (a) Synchronization region determined numerically: 
black represents parameters where < 10~^ for all the re- 
alizations of Xi{0), rjij and Tij. (b) Synchronization region 
determined from the stability analysis: black represents pa- 
rameters where the [Amaxl < 1 for all the realizations of rjij 
and Tij. The delays are Gaussian distributed with to = 3, 
— 100, all other parameters are as in Fig. 1. 



delays. 

Atay et al. have recently shown that with fixed 
delays the synchronization depends on the array archi- 
tecture: with the same total number of links, a random 
network exhibits better synchronization properties than 
a regular network. We investigated this issue in the case 
of distributed delays and found that the synchronization 
depends on the average number of links per node, < 6i >, 
but not on the architecture. We did simulations with ar- 
rays of different topologies TE] and found that the tran- 
sition to synchronization is similar for all the arrays pro- 
vided that they have the same < bi >. Figure 5(a) dis- 
plays the transition to synchronization as e increases for 
small- world and regular arrays with distributed delays, 
and for comparison. Fig. 5(b) displays results for the 
same arrays with fixed-delays. It can be observed that 
for distributed delays the transition to synchronization is 
independent of the array topology; however, it depends 
on the connectivity: the larger the number of mean links 
per node, the lower the coupling strength needed syn- 
chronize. In contrast, for fixed delays the synchroniza- 
tion depends not only on the connectivity but also on 
the architecture: for e large enough the arrays that have 
small-world topologies synchronize, but those that have 
regular topologies do not, in agreement with the results 

of a- 

Finally, let us assess the stability of the steady-state 
synchronized behavior found numerically by performing 
a statistical linear stability analysis. The delayed map 
Eq. can be written in non-delayed form by the in- 
troduction of a set of auxiliary variables 0j], yim(t) = 
Xi{t — m), where 1 < i < N and < m < M with 
Af = max(Tij). In terms of these new N{M -f 1) vari- 
ables Eq. Q becomes 
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FIG. 5: Influence of the architecture and the connectivity of 
the array for (a) Gaussian distributed delays (to = 3, c = 2) 
and (b) fixed delays {nj = 3 V i,j). Small-world topology 
with and an average of 10 links per node (o), 50 links per 
node (+); regular topology with 10 links per node (*); 50 
links per node (A). All other parameters are as in Fig. 1. 



yini{t + 1) = 

iyz,m-i{t) ifm^O, 
\ (1 - e)f[y,o{t)] + EU ^^^fbjj,r.. (*)] if ™ = 0. 

Next we define a vector of N{M + 1) components 
containing information about the present and past states 
of the array (for i = 1 to and for time ttot — M) : Y = 

(2/10,2/20, • ■ • ,2/Afo; 2/11, 2/21, ■ • ■ ,2/Afi; ■ • ■;2/iA/,2/2Af • ■ • ,2/JVAf) 
Eq. ^ can be re- written as Yi{t + 1) — 
Fi{Yi{t) . . .YN(^M+i){t)) and the synchronized state 
can be re- written as Y* — [xq . . . xq). 

Linearizing the equations of motion near Y* gives 
SYi{t) = Ej fijSYj{t), where the matrix F = (/y ) with 
fij = dFi/dYj can be cast as a set of {M + 1)^ blocks 
of dimension N x N. We denote these blocks as Tki, 
with k,l = 0, ...M. The blocks J^ki with fc > have all 
components equal to 0, except the blocks J-k+i,k which 
are N x N identity matrices. The elements of the blocks 
J^ki with A; = can be considered as composed by two 
parts: (i) (1 — e)f'{xo), which arise from the non-delayed 
term in Eq. and appear in the diagonal of ^00 ■ (ii) 
(e/&i)/'(a;o), which represent the contribution of a Hnk 
and due to the random delays appear in random posi- 
tions of the blocks. Each link contributes with two terms 
of the form (e/6i)/'(xo). There is a total of bi terms 
{e/bi)f'{xo) are which are distributed randomly in the 
blocks Toi- More precisely, they are located at the posi- 
tions and (j, i) in the blocks Ton^ and Torji respec- 
tively. The rest of the elements are zero. 

We calculated the eigenvalues of F for different real- 
izations of the connectivity and delay time distributions. 
The results are displayed in Fig. 4(b), where the black 
region indicates parameters where the maximum eigen- 
value of F, Xmax, has modulus less than 1 for all Tij and 
rjij realizations, and the white region indicates parame- 
ters where Xmax > 1 for all and rjij realizations. A 
very good agreement with the synchronization region de- 
termined numerically can be observed [the black region 



for c = observed in Fig. 4(a) does not appear in Fig. 
4(b) because in this region the synchronized dynamics is 
chaotic] . 

To conclude, we studied the dynamics of an ensemble 
of chaotic maps which interact with random delay times, 
and found that for adequate coupling strength the array 
synchronizes with all the maps in a steady-state. The 
synchronization is independent of the architecture of the 
array, but depends on the average number of links per 
node. We confirmed the numerical results by perform- 
ing a statistical linear stability analysis. Our findings 
provide another example of the nontrivial action of inho- 
mogeneities and disorder in coupled nonlinear systems; in 
our case, the presence of large enough disorder in the in- 
teraction times plays a constructive role, suppressing the 
chaotic dynamics of the individual units. We speculate 
that the " chaos-suppression by random delays" reported 
here might yield light into explaining the stable operation 
of many complex systems composed by nonlinear units 
which interact with each other with random communica- 
tion times. 
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